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Some time ago, Fisher and de Gennes pointed out that long-ranged cor- 
relations in a fluid close to its critical point T c cause distinct forces between 
immersed colloidal particles which can even lead to flocculation [C. R. Acad. 
Sc. Paris B 287, 207 (1978)]. Here we calculate such forces between pairs of 
spherical particles as function of both relevant thermodynamic variables, i.e., 
the reduced temperature t = (T — T c )/T c and the field h conjugate to the or- 
der parameter. This provides the basis for specific predictions concerning the 
phase behavior of a suspension of colloidal particles in a near-critical solvent. 
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I. INTRODUCTION 

In 1978 Fisher and de Gennes [0J predicted that the confinement of critical fluctuations 
of the order parameter in a binary liquid mixture near its critical demixing point T c gives rise 
to long-ranged forces between immersed plates or particles. In particular, they pointed out 
that these long-ranged forces would eventually lead to the flocculation of colloidal particles 
which are dissolved in a near-critical binary liquid mixture 0. A few years later, such a 
solvent-mediated flocculation was observed experimentally for silica spheres immersed in the 
binary liquid mixture of water and 2,6-lutidine H,|J. However, the precise interpretation of 
these experiments is still under debate; although fluctuation- induced forces as predicted by 
Fisher and de Gennes certainly play a major role, additional mechanisms, such as screening 
effects in the case of charged colloidal particles 0, may also contribute ||. In the meantime 
additional experimental evidence for this kind of flocculation phenomena has emerged for 
other binary mixtures acting as solvents . 

From a colloid physics perspective effective, solvent-mediated interactions between dis- 
solved colloidal particles are of basic importance p|-[12|. The richness of the physical prop- 
erties of these systems is mainly based on the possibility to tune these effective interactions 
over wide ranges of strength and form of the interaction potential. Traditionally, this tun- 
ing is accomplished by changing the chemical composition of the solvent by adding salt, 
polymers, or other components []9|— |T2 1 . Compared with such modifications, changes of the 
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temperature or pressure typically result only in minor changes of the effective interactions. 
On the other hand, however, effective interactions generated by bringing solvents close to a 
phase transition of their own are extremely sensitive to such changes. In this broader con- 
text, flocculation of colloidal particles induced by critical fluctuations near a second- order 
phase transition of the solvent as predicted by Fisher and de Gennes ]T]|| provided the 
first instance in which effective interactions are tuned by changing the physical conditions 
rather than the chemical composition of the solvent. Later on, this idea has been extended 
to first-order transitions of the solvent, for which wetting transitions can occur which give 
rise to wetting films of the preferred phase coating the colloidal particles |T^TB|. In the 
meantime the study of effective interactions between dissolved particles generated by chang- 
ing the physical conditions of the solvent, and the resulting phase behavior of the colloidal 
suspension, has emerged as one of the main fields of research in colloid physics [13 . 

The physical origin of the force predicted by Fisher and de Gennes is analogous to the 
Casimir force between conducting plates, which arises due to the confinement of quantum 
mechanical vacuum fluctuations of the electromagnetic field |T7|-[l9|. (Related phenomena 
occur in string theory 



liquid crystals (see, e.g., ref. ||21||), and microemulsions |22 
However, in contrast to the present so-called "critical Casimir effect" , in liquid crystals the 
director-fluctuation induced forces are notoriously difficult to separate from the background 
of dispersion forces due to the absence of a singular temperature dependence P^| .) In recent 
years this "critical Casimir effect" has attracted increasing theoretical |2^-|28| and exper- 
imental P J29| , |30[| interest. In this case the fluctuations are provided by the critical order 
parameter fluctuations of a bona fide second-order phase transition in the corresponding 
bulk sample. For He 4 wetting films close to the superfluid phase transition the theoreti- 
cal predictions |25 for critical Casimir forces between parallel surfaces exhibiting Dirichlet 
boundary conditions have been confirmed quantitatively JHJ. These quantum fluids offer 
the opportunity to study critical Casimir forces in the absence of surface fields, i.e., the 
purely fluctuation induced forces. 

In a classical binary liquid mixture near its critical demixing point, the order parameter 
is a suitable concentration difference between the two species forming the liquid. In this 
case the inevitable preference of confining boundaries for one of the two species results in 
the presence of effective surface fields leading to nonvanishing order parameter profiles even 
at T > T c [PH| . (This is the case discussed by Fisher and de Gennes These so-called 

"critical adsorption" profiles become particularly long-ranged due to the correlation effects 
induced by the critical fluctuations of the order parameter of the solvent. In the case of 
a planar wall, critical adsorption has been studied in much detail fl32|-j39|j . Asymptotically 
close to the critical point T c it is characterized by a surface field which is infinitely large so 
that the order parameter profile actually diverges upon approaching the wall up to atomic 
distances. Later on, critical adsorption on curved surfaces of single spherical and rodlike 
colloidal particles has been studied, where it turns out, in particular, that critical adsorption 
on a microscopically thin "needle" represents a distinct universality class of its own |f40| , f4ll . 
Critical adsorption on the rough interface between the critical fluid and its noncritical vapor 
occurring at the critical end point of the binary liquid mixture has been addressed as well 



|3| ; likewise the effect of quenched surface modulations and roughness [ |3 . The interference 



of critical adsorption on neighboring spheres gives rise to the critical Casimir forces which 
have been argued to contribute to the occurrence of flocculation near T c P,fH|-r4S|] . However, 
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a quantitative understanding of these phenomena requires the knowledge of the critical 
adsorption profiles near the colloidal particles and the resulting effective free energy of 
interaction in the whole vicinity of the critical point, i.e., as function of both the reduced 
temperature t = (T — T c )/T c and the bulk field h conjugate to the order parameter. So 
far, this ambitious goal has not yet been accomplished. Instead, the introduction of the 
additional complication of a surface curvature has limited the knowledge of the corresponding 
critical adsorption so far to the case of spheres for the particular thermodynamic state 
(t, h) = (0, 0) of the solvent @g7|-||. 



In a previous work [iB|, at least the temperature dependence of the critical Casimir force 
between a sphere and a planar container wall has been discussed. It turns out that the 
force becomes maximal not at but above T c . In ref. [|(| also the force F between a pair of 
spheres of radius R at distance D + 2 R between their centers (see Fig. [I]) has been briefly 
analyzed. Close to the critical demixing point, this force separates into a regular background 
contribution and a singular contribution -F s f ng of universal character, which is attractive if 
the same of the two coexisting bulk phases is enriched near both spheres. The results for 
F sing can be cast into the form 



F± ng (D } R;t,h) = J^ ± (A = -,e ± = -,H = sgn(/i)- 

where £ = £± is the true correlation length governing the exponential decay of the order 
parameter correlation function in the bulk for t = (T — T c )/T c ^ at h = 0; Ih is the cor- 
relation length governing the exponential decay of the order parameter correlation function 
for h — > at t = 0. K± are universal scaling functions. Note that surface fields hi break 
the symmetry h — > —h, i.e., F s f ng depends on the sign of H. Here and in the following we 
consider only the fixed point value hi — > +oo. 

The aim of this work is to study the critical Casimir force between two spherical colloidal 
particles as function of both relevant thermodynamic variables t and h. Figure || shows the 
schematic phase diagram of the solvent, e.g., a binary liquid mixture. Different thermody- 
namic paths are indicated along which the dependence of the force on t and h is studied. 

The remainder of this work is organized as follows: In Sec.|TJ we define the field-theoretic 
model used to describe the near-critical solvent. In Sec.|T| the numerical calculation of 
the adsorption profiles in mean-field approximation is discussed. In Sec. we present the 
results for the force between two spherical particles and compare them with the limiting 
cases of small and large particle radii. The discussion of the results for two parallel plates, 
which inter alia are relevant for a pair of spheres via the Derjaguin approximation, is left to 
Sec. [V[ Finally, in Sec. VI, we summarize and discuss our results with regards to the phase 



behavior of a suspension of colloidal particles in a near-critical fluid. 



II. FIELD-THEORETIC MODEL 



Near criticality the behavior of the system is governed by fluctuations of the order param- 
eter $ on a large length scale such that the emerging universal properties are independent 
of microscopic details. A coarse-graining procedure leads to the continuum description in 
terms of the Ginzburg-Landau Hamiltonian [p^, p3| , p0[| 
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W{*(r)} = J d<V{l[V$(r)] 2 + ^$ 2 (r) + ^ $ 4 (r) - h*(r)} (2.1) 

v 

of an order parameter field $(r) in the bulk of volume V. For the liquids considered here 
$(r) is a scalar field. The Hamiltonian determines the statistical weight exp[— ^{^(r)}] of 
the configuration {$(r)}. The external bulk field h(r) is conjugate to the order parameter 
and breaks its symmetry. The temperature dependence enters via 

T ~t = —^ . (2.2) 

-L c 

The partition function of the system is given by 

2(r l / l ;«)=/p{$}exp[-W{$(r)}] , (2.3) 



where denotes an appropiately defined functional integration |33|. The corresponding 

singular part of the free energy reads 

•^sing^, h; u) = -k B T c In Z(t, h; u) . (2.4) 

The coarse-graining procedure leading to a continuum model from a given lattice model 



can be generalized to systems with surfaces resulting in the Hamiltonian [p2| , p3 



nmr)} = J ^r{i[V$(r)] 2 + ^$ 2 (r) + ^$ 4 (r) - fe*(r)} 
v 

+ J dS{U 2 (v s ) - ht^Tg)] . (2.5) 

s 

The first integral runs over the volume V available for the critical medium, which is bounded 
by the surface(s) S — dV. Apart from the restricted space integration, the first integral in 
Eq. ( |2.5|) has the same form as the Hamiltonian fl2.1|) for bulk systems disregarding surfaces. 
The second integral in Eq. ( |2.5| ) represents a surface contribution in which the integral runs 
over the boundary surface S = d V; c is related to the strength of the coupling of critical 
degrees of freedom near the surface and hi is the surface analogue of the bulk field h |3"2"| , |3"3" 



The systematic treatment of the critical fluctuations as implied by Eq. ( |2.3| ) and a jus- 
tification of the form of 7Y{$(r)} is provided by the field-theoretic renormalization group 



framework [33]. The basic idea is to consider the physical system on larger and larger length 
scales and to single out those quantities which eventually remain unchanged under the re- 
peated action of such scale transformations. Such quantities tend to their fixed point values 
right at the critical point |50[| . For bulk systems, there remain only two relevant adjustable 
parameters, corresponding to r and h. The surface generates a subdivision of a given bulk 
universality class into different surface universality classes characterized by different fixed 
point values for c and hi. For hi = the corresponding leading singular behaviors are 
determined by the fixed-point values c = +oo for the ordinary transition (O), c = for the 
special or surface-bulk transition, and c = — oo for the extraordinary or normal transition 
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(E) []32 , |33 . For the O and E transitions only r, h, and hi remain as relevant variables. The 



three cases above were originally studied for magnetic systems. For liquids, however, E is 
the generic case and therefore it is called the normal transition. The case hi = oo with c 
arbitrary corresponds to critical adsorption, but the leading singular behavior is the same 



as for the E transition 36 



Neglecting fluctuations the Ginzburg-Landau Hamiltonian 7i{$(r)} itself represents a 
free energy functional which upon minimization with respect to $(r) yields the Euler- 
Lagrange equations for order parameter profiles and correlation functions. The minimization 
is equivalent to the mean- field description of the system, 



6H{®(r)} 



= , (2.6) 



where the fluctuating order parameter $ is replaced by its mean value ($). The universal 
critical behavior is captured correctly by mean-field theory if the space dimension d exceeds 
the upper critical dimension d uc = 4. 

Because of the relative complexity of the geometry considered here (see Fig. [I]), it is 
not possible to carry out the complete renormalization procedure in closed form. Instead, 
available knowledge for bulk systems in d = 3 is used in combination with the mean-field 
solution for the geometry considered here in order to estimate universal quantities in d = 3. 
An alternative way to overcome mean-field theory is followed within the concept of the 
Local Functional Theory developed in refs. |5l|^4|, in which the free energy is described by 
a functional depending locally on the order parameter and its derivative. With mean-field 
theory as a special case, it takes into account available knowledge in d = 3, such as the 
correct values for the critical exponents. 



III. ORDER PARAMETER PROFILES 

In this Section we describe some technical aspects of the tools and numerical methods 
used for obtaining the mean-field results, i.e., the determination of the order parameter 
profiles at the critical point, the short distance expansion of the profile close to a surface, 
the minimization procedure to determine the profiles away from the critical point, and the 
calculation of the force between the particles based on the stress tensor. 

Minimization in Eq. ( |2.6| ) leads to the differential equation 

- Am + rm + m 3 - H = , (3.1) 

where the coupling constant u is absorbed in the order parameter and in the bulk magnetic 
field, 

m=y|<*>,ff=y|/i. (3.2) 

Note that we are interested in an estimate of the order parameter profile m in d = 3 where, 
at the fixed point, u is a positive number. For d /* 4 the mean- field solution is exact, but 
the fixed point value of the coupling constant u vanishes in this limit. 
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The surface integrals in Eq. ( |2.5|) lead to a boundary condition for the order parameter 
at the surfaces of the spherical particles. At the critical adsorption fixed point this boundary 
condition is given by 

$ s = oo . (3.3) 



A. Order parameter profiles at criticality 

While our main interest is the dependence of the Casimir force on t and h, the knowledge 
of the order parameter profile at the critical point, i.e., for (t, h) = (0, 0), is a useful starting 
point. In ref. this latter profile is given for a critical system confined to a spherical shell 
between two concentric spheres with radii R± at which the order parameter diverges. In 
terms of elliptic functions it reads 




m(r) = v — — — , (3.4) 

r cn 




where r is the radial distance from the common center of the spheres and k is the module of 



the elliptic function cn [55 1. The module k must be adjusted such that the profile diverges 



at the given radii R± leading to the following implicit equation for k: 



^=(f t ) 2 = exp(±v/2P^lA^)), (3.5) 



where R = ^/R+ R- is the geometric mean of the radii and K is the complete elliptic 
integral of the first kind. 

In ref. MM it is outlined how the results for concentric spheres can be used for different 



geometries via a conformal transformation of the coordinates, 

r' r + R R 

~ |r + R| 2 ~ 2IRI 2 ' 



(3.6) 



where r' and r are the new and the old coordinates, respectively, and R defines the symmetry 
axis and the eventual shape of the new geometry. For |R| = i? the two concentric spheres 
with radii R± are mapped onto two spheres of equal size which is the generic case we are 
interested in. For |R| = R + the two concentric spheres are mapped onto a single spherical 
particle in front of a planar wall for which results were obtained in ref. MM . 

Under the conformal transformation (|3.6| ), each local scaling field is multiplied by a scale 
factor 

R r' Ir'l 2 
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which is raised to the power of the scaling dimension of this scaling field. In mean-field 
approximation and d = 4, the scaling dimension of the order parameter $ is 1, which 
implies for the order parameter profile 



m r 



" -Kr'r^coBctKr')!), (3.8) 



where ?^ conc (|r|) is the profile in the concentric geometry given by Eq. 

Let D denote the surface-to-surface distance between the two spheres and R their com- 
mon radius; then the corresponding radii R± in the concentric geometry are given by 



D . ( :\ ■ 1 X =M 



^ = T A IJ3t) (3-9) 



with the abbreviation 



A = : 1 • 1^. (3.10) 



Putting the pieces together, Eqs. (|3.2|) - (|3.9|) yield the order parameter profile in the geometry 
of two spherical particles with radius R a distance D apart from each other right at the 
critical point (t,h) = (0,0). 



B. Short distance expansion 



In the following we describe the methods for obtaining results off criticality, i.e., for 
t ^ and h ^ 0, i.e., r ^ and H ^ 0. Equation ( jSTlf ) will then yield a different profile 
m(r'). However, the profile at the critical point serves as a suitable starting point for the 
minimization procedure (2J5). Close to the surfaces it is even possible to calculate the 
deviation of the order parameter from the profile at the critical point analytically. In this 
context the so-called short distance expansion turns out to be useful, and has been carried 

With H 7^ in addition, it reads: 



out for r ^ in ref. |40 



m(s -> 0;R, t,H) 



V2 V2d- 1 



6 R 



(3.11) 



y/2 
6 



S< — T 



5 {d 



IQ R? 



(d- l)(d-2) 
R? 



y/2 2 / H d-l Ad 3 - 39d 2 + 120d - 85 
+ ~T S 1^2 _ SR T 54 R 3 

+ 0(s 3 ) , 



where s is the radial distance from the surface of a sphere with radius R. The difference 
5m = m(s — ► 0; R, r, H) — m(s — ► 0; R, r — 0, H — 0) gives the deviation from the profile at 
the critical point close to the surfaces. While the profile itself diverges at the surfaces, the 
deviation 5m vanishes with the radial distance s as 
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\/2 

5m(s — ► 0; R, r, if) = r s 

6 




v/2d- 1 
~3 R~ 



s 2 + 



(3.12) 



where the bulk field H and the radius R of the spherical particles appear only in 0(s 2 ) 
and higher order. Equation (|3.11| ) is only valid for small s since the presence of the second 
particle is neglected. Corrections to Eqs. ( |3.11| ) and (3.12), however, are of order 0(s 3 



as 



derived for the case of parallel plates, where they are called distant wall corrections |H,|56|. 



C. Numerical solution for the complete profiles off criticality 

So far we have provided results for the order parameter profile at the critical point for the 
two sphere geometry [Eq. ( ^.8|) 1 and off criticality but close to the surface of a single sphere 
[Eq. (|3.11|) 1- In order to obtain the complete order parameter profile away from the critical 
point for the geometry consisting of two spheres with equal size (see Fig. |1]) additional, 
numerical effort is required. 

Upon discretization the order parameter profile m(r'; D, R; r, H) in the two sphere ge- 
ometry (Fig. [1]) is described by a set of parameter values where i and j indicate the 
dependence on the two independent spatial variables required for the present cylindrical 
symmetry around the axis between the two centers of the spheres [see, c.f., Eq. ( p.l5| )j: 

m(r';D,R;r,H) = m{{a ij }) . (3.13) 

Within mean-field theory Ti.{{a,ij}) has to be minimized with respect to {a^} for given D, 



R, r, and H . To this end we apply the method of steepest descent f57[j 



(n+l) (n) dH({aij}) 



(n) 

a- 



(3.14) 



» 



where k is a positive number and is the n-th iteration of the parameter a^. 

It is known that in general the method of steepest descent converges slowly. Nevertheless 
in the present case it is found to be useful since the dependence on the choice of the initial 
set of parameters {a^f 1 } is weak. In fact, it turns out that the profile at the critical point is a 
rather useful starting profile. The convergence also depends on the value of k: if it is chosen 
too large, the convergence breaks down, if it is chosen too small, the convergence becomes 
too slow. As we will discuss later, one advantage of the present minimization method is 
that the explicit evaluation of the Hamiltonian 7i is not necessary, since only the gradient 
d^ii"'^ ^ s nee ded, which can be obtained independently. In order to keep this advantage, 
iterative methods to determine the best choice for the value of k, which explicitly need the 
value of the functional that is to be minimized, are not applied. Nonetheless the best choice 
of k is the largest value for which the minimization still proceeds. 

In particular, we use the parametrization 

m(4; D, R; r, H) = m(r' ij - 1 D, R; 0, 0) + a {j , (3.15) 

where y[- denotes discretized spatial positions in the two sphere geometry, ^(r^; D, R; 0, 0) 
is the value of the order parameter profile for (r, H) = (0, 0) at those positions, according 
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to Eq. ( |3.8|) , and denotes the deviation from this value at the position used in the 
minimization procedure ( 3.14Q . Thus starting the procedure with the profile for {t,H) = 
(0, 0) implies the choice 

4 0) = 0, Vi,j. (3.16) 



In addition, Eq. ( |3.12| ) provides the behavior of the for r^- close to the surface, i.e., 



aff = for all n if r/. = r^. From Eq. (|3.12|) even the slope of a function f(s) used for 



interpolating the as function of the radial distance s from the surface of the sphere is 
known for s — > 0. 

The cylindrical symmetry of the problem allows one to consider the planar geometry 
depicted in Fig. 0. Instead of discretizing the space by a square lattice we use a lattice which 
takes into account that there with strong changes in the order parameter m, i.e., 

close to the surfaces, and areas with weak changes, i.e., far away from the spheres. In those 
areas with weak changes less lattice points are needed. For this purpose we use what we call 
the conformal lattice. The idea is to use the conformal mapping ( ^.6[ ) in order to create this 
lattice from that one for the case of concentric spheres. In analogy to electrostatics the field 
lines and the lines of equal potential of a system with two concentric spheres forming the 
two electrodes of a capacitor are mapped onto the two sphere geometry such that the new 
lines are the field lines and the lines of equal potential of a system with two charged spheres. 
The intersections of these lines form the lattice sites of the conformal lattice. The advantage 
of this lattice is that the field lines approach the spheres perpendicularly to their surface, 
which allows one to use the short distance behavior [Eqs. ( p. 11 ) and ( |3.12|) 1 as a boundary 



condition for the a^ . Figure [3] shows the conformal lattice for the geometry of two spheres. 
In principle, it would be possible to use only one half of this lattice, since this system is 
symmetric with respect to the vertical midplane. However, this reduction leaves one with 
a disadvantage. The value of the order parameter at this midplane wall is not known, and 
thus it must be determined within the minimization process. The minimization turns out 
to be inefficient then. Therefore for the two sphere geometry we resort to the full conformal 
lattice enclosing both spheres as shown in Fig.|3|. 
The Hamiltonian 7i [Eq. ( |2.5|) 1 is given by 

U = - [ dVCim) + - I dSCsim) , (3.17) 
u J u J 

V s 

with 

1 r 1 

C(m) = -(Vm) 2 + -m 2 + -m 4 - Hm (3.18) 

and Cs(m) = ^m 2 — Him with Hi = hi Performing the functional derivative with 

respect to the order parameter m leads, via partial integration, to 



5H{m(r)} 6 



5m(r) 



- [ dV{-Am + Tm + m 3 - H} + - I dS {-d n m + C' s (m)} , (3.19) 
u Jv u J s 



where d n m is the derivative of the order parameter in direction of the surface normal. We 
note that the surface integral vanishes, if the boundary conditions [Eqs. ( |3.11| ) and ( |3.12|) 1 
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are fulfilled. Equation ( |3.14|) contains the derivative of the Hamiltonian Tt with respect to 
the parameters which represent the order parameter m. With m = m(aij) = const in a 
small volume 8V{t^) around the position the gradient in Eq. ( |3.14j ) is approximated by 

dH ^ ai ^ - 6 8V(r'){-Am + Tm + m 3 -H} , , . (3.20) 



d a,ij u 



r = r 



Note that a constant prefactor in Eq. ( |3.20| ) can be absorbed by a rescaling of k in Eq. ( |3.14j ). 
If Eq. (|3.1|) holds, i.e., if the expression in curly brackets in Eq. ( 3.20 ) vanishes, Eq. ( |3.11| ) 

reduces to 

at 1] = 4 B) > ( 3 - 21 ) 

i.e., the minimization has led to a fixed point solution of the iteration. 

For the minimization process the second derivatives of the order parameter profile are 
needed [see Eq. ( |3.20| )1. Since the order parameter is calculated only at the discrete sites 
r-j, numerical differentiation must be applied. For a numerical derivative, in general, it is 
preferable to have the sites on an equidistant square lattice. However, it is also possible 
to deal with derivatives on the conformal lattice. To this end we have used the spline 
approximation |58| in order to interpolate between all the on each of the lines of the 
conformal lattice (see Fig. |3|). This approximation directly yields the second derivative of 
the order parameter at the lattice sites with respect to the variables parametrizing the lines. 
Under a conformal mapping angles are conserved, so that the lines in the conformal lattice 
intersect perpendicularly which allows one to translate the derivatives obtained from the 
spline approximation into derivatives with respect to the cylindrical coordinates r' = (p, z). 



D. Stress tensor 

With the tools now available to determine the order parameter profile in the two sphere 
geometry the next step is to calculate the effective forces between the two spheres. The 
straightforward way to determine such forces is to calculate the free energy —ksT \nZ as 
function of the distance D between the spheres, and then to estimate the derivative with 
respect to the distance by considering finite differences. However, the total value of the 
free energy is expected to be large as compared with the differences, and thus one has to 
deal numerically with differences between large numbers, which causes high inaccuracies. 
Therefore we resort to an alternative method, in which the force is calculated directly using 
the so-called stress tensor Tu U , which is defined by the linear response 



6H 



d d r ^ T, 



» 



(3.22) 



where b is a non-conformal coordinate transformation, \i and v index the spatial coordinates, 
and 57i is the corresponding energy shift. With the Lagrangian C(m) in Eq. QTT8 ) it reads 



44,46.5 



T, 



6 

u 



o{o u m) 



(3.23) 



10 



The knowledge of the stress tensor on any infinitely extended surface which separates the 
two spheres from one another, or likewise on any closed, finite surface which contains only 
one of the spheres, allows one to calculate the force between the spheres. This is illustrated 
in terms of the spherical volume Vs with radius R containing the right sphere in Fig.|j. 
Using the coordinate transformation 



b z = a, Vr e V s 
= 0, Vr £V S 



(3.24) 



the volume V$ is shifted in z-direction by the amount a. The derivative ^ L then leads to 
a delta function such that the integral over the whole volume V reduces to an integral over 
the boundary of Vs- In this case the force F s - mg = kB ^ c K with the universal scaling function 
K [see Eq. is given by 



dm 

da 



da J v d r 



ldV s 



dip sin 2 <p [T ZZ (R, <p) cos <p + T zp (R, cp) sin ip] 



(3.25) 



where n is the unit vector perpendicular to the boundary dVs- The quantities used in the 
parametrization in Eq. ( |3.25| ) are sketched in Fig. |j. 

In the geometry with cylindrical symmetry [m = m(p, z)\ considered here, Eq. ( ^.23| ) 
leads to 



T 

1 7'. 



1 

u 



/ dm 
^\dl 



+ 4 



, m dm 
P dp 



+ 2m- 



d 2 m 



dm 
dp 



3rm A 



-m 



(3.26) 



and 



T. 



zp 



1 



I dm \ / dm \ 
dz ) \dp J 



2m 



d 2 m 
dpdz 



(3.27) 
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where the so-called improvement term — ^(V^V^m 2 — 5 pu Am 2 ) has been added 

In our calculation, we choose Vs as indicated in Fig. |3], determine the stress tensor at 
the lattice points, and use spline interpolations between them. Note that the possibility of 
choosing different spherical surfaces enclosing one of the spheres, which via Eq. ( 3.25| ) should 



all result in the same force -F sing , provides a valuable and sensitive check for both the validity 
and the accuracy of the numerical results. 



IV. RESULTS FOR THE TWO-PARTICLE INTERACTION 



The goal of this work is to provide a major step towards the understanding of aggregation 
of colloidal particles dissolved in a binary liquid mixture close to criticality. In the previous 



11 



section we explained the various methods for obtaining the critical Casimir force between 
two spherical particles of equal size. The results are presented in this section. 

The aggregation phenomenon is driven by the dependences on the temperature t = 
(T — T c )/T c and on the bulk composition ca of the solvent. While the former is related to 
the temperature parameter r the latter is related to the bulk field H . Thus also the Casimir 
force between two particles must be determined as a function of those two parameters. We 
present results for the scaling function K±(A = ^, B± = H = sgn(h) j^) which is related 
to the force F s f ng between the particles via Eq. fll.lj) . The correlation length £± for H = 
depends on the temperature T as [see also Eq. ( |2.2|) 1 



with v — | within mean-field theory while the correlation length In at T c depends on the 
bulk field H as 



— 



ca - ca,, 



-v/P 



1 



-^|F|- 1/3 , r = 0, (4.2) 



CA,c 

where u//3 — 1 within mean-field theory; £ and Zo are nonuniversal amplitudes. The first 
parts of Eqs. (|4.1| ) and fl4.2|) hold in general whereas the second parts correspond to the 
present mean-field description. The scaling variable 5 = sgn(H)D/lH depends on the sign 
of H [Eq. (tL.l|) 1. The relation between the bulk field H and the concentration ca will be 
discussed further in Subsec. [IV GL but we keep in mind that a negative value of H means that 
the fluid component, which is preferred by the particles, is the one with a low concentration 
in the bulk (see Fig. 0) . 



A. Dependence on temperature at the critical composition 



The variation of the temperature leads to an interesting behavior of the force [46 . 
According to Fig.|5] the force exhibits a maximum for a temperature above the critical 
point. The position of the maximum of the force as function of the scaling variable @ + 
depends on the distance D between the particles. In the limit of large distances the position 
©max = @max(A, H = 0) of the maximum approaches that one of the critical point, i.e., 

9 max (A = ^ ooj . (4.3) 

On the other hand, for small distances the position of the maximum approaches a finite 
value above T c 

e max ( a - o) - eL°L . (4.4) 

These two limiting behaviors are shown in Fig. |5|. The disappearance of the maximum for 
large distances along with the derivation of the small radius expansion will be discussed in 
more detail later on in this section. 
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B. Dependence on composition at the critical temperature 

The above results show that there are rich structures in the force curves off criticality. 
Therefore it is important to monitor the force not only at the critical point, but also in its 
vicinity. This becomes even more evident if not the temperature but the composition of 
the solvent is varied. Along the dashed path in Fig. ^| and for small distances D the force 
exhibits a pronounced maximum as can be seen in Fig. [| For large distances there is only 
a maximum at the critical point, i.e., H = 0. For H > the force decays monotonously for 
all distances D. 

Figure [j] shows the maximum values of the force for a given distance D as function of the 
temperature (solid line in Fig.|2|) or the composition (dashed line in Fig.|2[). For comparison, 
also the force at the critical point for the same particle distance D is plotted (dash-dotted 
line in Fig.|?p. For intermediate and larger distances the three lines merge, which is due to 
the fact that for large distances the position of the maximum of the force approaches that 
of the critical point. For small distances the differences are pronounced. Obviously, the 
composition of the solvent plays an important role. 



C. Dependence on temperature off the critical composition 

So far the thermodynamic paths were chosen such that either the temperature or the 
solvent composition were fixed at their critical point value. However, it is experimentally 
relevant to consider also thermodynamic paths along which the temperature is varied at fixed 
compositions ca ^ ca,c- The corresponding results are given in Fig. || for a fixed distance 
D. For small deviations \ca — ca,c\ ^ °a,c the position of the maximum of the Casimir force 
is more or less unchanged while its absolute value is changed considerably. However, the 
main features of the temperature scan do not differ much from that at ca, c so that they are 
sufficiently robust in view of possible experimental tests with finite resolution for ca- 

If the dotted path in Fig. |2| is extended to temperatures below T c the force is analytic as 
function of the temperature until the path intersects with the coexistence line. 



D. Derjaguin approximation 

In the limit that the radius R of the spherical particles is much larger than both the 
correlation length £ and the distance D of closest approach surface-to-surface between the 
two particles, the particles can be regarded as composed of a pile of fringes. Each fringe 
builds a fringe-like slit with distance L = D + rjj/ R, where r\\ is the radius of the fringe. 
Within this Derjaguin approximation and in d = 4, the force between the particles in 



units of [see Eq. is given by 

■5/9 f°° 9 ^±(6±(l + i; 2 ),S(l + i; 2 )) 

K ± A,6±,S =4vrA- 5 / 2 / dv v ±K ±V h K >± 

Jo 



(4.5) 



where v is a dimensionless integration variable and K±(L/£±,sgn(H)L/lH) is the analogous 
scaling function for the Casimir force in the slit geometry with parallel walls at distance L 
[see, c.f., Sec.0. 
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E. Small radius expansion 



If on the other hand the radius R — albeit large on the microscopic scale — is much 
smaller than D and £, the statistical Boltzmann weight e~ SHs characterizing the presence 
of the sphere centerd at r$ can be systematically expanded in terms of increasing powers of 
R i.e., 



e - SHs /d^ = l + tfR x * $(r s ) + cf iT* 2 $ 2 (r 5 ) + . . . (4.6) 

where is an amplitude and where x$ = f3/v and x$2 = d — v~ l (with the standard 
bulk critical exponents (3 and v) are the scaling dimensions of \1/ = $,$ 2 . The ellipses 
stand for contributions which vanish more rapidly for R — > 0. The coefficients are fixed 
by c* = A^/By, where and By, respectively, are amplitudes of the half-space (hs) 
profile {^(z))l s T=Tc = [2z)~ x ^ at the critical point of the fluid for the boundary condi- 
tion | corresponding to critical adsorption, and of the bulk two-point correlation function 
(^(r)^(0))t,,T=T c = B^r~ 2xii P4] |. Since both spheres are equally small the free energy T 



follows from the bulk average of two statistical weights of the type given in Eq. 

> e -SH Sle -6Hs 2 } h 



r{D)-n<*>) _ ln 

k B T c 



(4.7) 



where the spheres S% and S2 are a distance D apart from each other. Using Eq. ( |4.6|) we 

find [||; 



F{D) J^(oo) _ , $ pxft^/flj,,, . , ,\, n _ .-,,.<i> /? .-, , 
knT c 1 



-(4 R<°*)*($(r Sl ) $(r 52 )) b [1 - 2c* + ...], (4.8) 



where ($(r5 1 ) $(rs 2 ))b is the bulk two point correlation function and ($)b provides a measure 
of the deviation of the bulk composition from its value at the critical point. The Casimir 
force is the derivative of the free energy with respect to the distance D between the spheres. 



Note that in Eq. (|4.8|) the first term in the square brackets gives rise to an attractive force 
which is increased by the second term if ($) b is negative, i.e., if the binary liquid mixture 
is poor in the component preferred by the colloids (compare Subsecs. |IVB| and [IV Q ). 

Within mean-field approximation, the correlation function is given by (see Eq. (B5) in 
ref. 0) 

($(r 5l )$(r S2 )) 6 = G+(|r 1 -r 2 |= J D;C + ) 



1 2 2 - d ' 2 ( D 



B *T^~ 2 ( - ) K^(D/^ + ) , (4.9) 



1 , T -> T, 



c 



where K(d_2)/2 is a modified Bessel function and where £ + = Ih for the case r = and 
H 7^ 0. Thus as expected two point-like perturbations of the bulk do not exhibit a rich 
interaction structure but a monotonous decay. If a small sphere is immersed at a distance 
D from another small sphere with a given adsorption profile, the resulting perturbation of 
this profile leads to an effective interaction. At those distances, at which the small radius 
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expansion can be applied the profile of the distant sphere is already close to the bulk value. 
This differs from the system studied in ref. [46] where a small sphere in front of a planar 



wall is considered. In that case the sphere is exposed to the half-space adsorption profile 
which is a much stronger perturbation of the bulk order parameter than that of a sphere 
(compare refs. [HO] and pl|j). 



F. Comparison with dispersion forces 

It has been stressed above that the Casimir forces depend sensitively on both temperature 
and composition of the solvent. The Casimir forces add to the omnipresent dispersion 
forces, which also depend on the thermodynamic state of the solvent. However, this latter 
dependence is smooth and - within the critical region of interest here - it is weak. Therefore 
the dependence of the Casimir force on temperature and concentration allows one to extract 
it from the actual total force which is the only one which is experimentally accessible. As 
an example for a binary liquid mixture we consider a water-lutidine mixture for which the 
critical temperature is T c = 307.1 K ||. The strength of the dispersion forces is characterized 
by the Hamaker constant, which for a system of silica spheres immersed in such a mixture is 
approximately A = 10 -20 J. Figure^ shows the comparison of the dispersion forces and of 
the Casimir forces at the critical point for spheres with a typical radius R = 100 nm which 
add up to the total force 

F tot (A, 0, 0) = K+(A, 0, 0) + F disp (A) . (4.10) 



The dispersion forces are given by 



A 32 1 
F disp (A) = "^y ( A + 2 ) ( A2 + 4 A + 4 ) ( A + 4 )2 A 2 ' ( 4 - n ) 

where A is the aforementioned Hamaker constant. In the limit of small distances the Casimir 
force is known from the Derjaguin approximation 

K+(A -> 0,0,0) = -vrA n A~ 2 , (4.12) 



where A-^ = 0.326 is the estimate of the Casimir amplitude in d = 3 [£7j], and in the limit 



of large distances the Casimir force can be obtained from the small radius expansion [4T5| 



K + (A oo,0,0) = -2^ TT (A + 2)- 2/3 /^ 1 , (4.13) 



where (5/v = 0.52 and A n = = 7.73 || 



In Fig. |S] also the forces between a single sphere and a planar wall are shown, as considered 
in ref. [|46[]. The dispersion forces then read F disp (A) = — ^| ( A+2 1 - ) 2 A 2 , where A is the surface- 
to-surface distance in units of the radius R. For small distances the amplitude of the Casimir 
force is two times bigger than in the case of two spheres while for large distances the small 



radius expansion leads to [49 



K + {A -> oo,0,0) = -2^ TT (2A)- /3/ "- 1 . (4.14) 
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For small distances all forces increase as A -2 , only the amplitudes differ from each other. 
The values considered for A and T c lead to Casimir forces that are about ten times larger 
than the dispersion forces. For large distances, however, the (unretarded) dispersion forces 
(~ A~ 7 ) decay faster than the Casimir forces (~ A -2 04 ). (Ultimately retardation of the 
dispersion forces leads to a decay ~ A~ 8 and ~ A~ 5 for SS and SW, respectively. For A = 1 
the Casimir forces between the two spheres are about fifty times larger than the dispersion 
forces such that the latter contribution to the total force [Eq. ( [4.10 )] can be neglected. For a 
sphere in front of a planar wall and for large distances the Casimir forces (~ A~ L52 ) decay 
even slower than they increase for small distances (~ A -2 ). 



G. Equation of state 



For the demixing transition of a binary liquid mixture the difference {^a—^b) — {.^a—^b)c 
of the chemical potentials of the two components is proportional to the bulk field h. The 
potential difference leads to a concentration ca 7^ ca, c of species A in the bulk. The bulk 
order parameter Mb is proportional to the concentration difference Mb = A(ca — ca, c ), where 
the amplitude A depends on the actual definition of the order parameter. For such a system 
in the vicinity of its critical point the equation of state is given by [^3 



h = VM 5 b f\t 



B 



M h 



1//3 N 



(4.15) 



where j3 and 5 are the standard critical exponents, T> and B are nonuniversal amplitudes, 
and the scaling function f(x) — 1 + x describes the crossover between the critical behavior 
at t = and h = 0, respectively. Eq. (|4~T5T) yields M b = ±B\tf for h -> and M b = (^f/ 5 
for t — > as limiting cases. The choice of the amplitude T> defines a field h. Since, however, 
the combination ^ dVhMb is the field contribution of the statistical weight of the system, V 
is fixed by the amplitude of the chosen order parameter. 

Away from the critical point the two point correlation function decays for large distances 
r exponentially as Gb(r) ~ exp(— r/^)/r^ -2 ^ 2 , with the true correlation length £. Within 
mean-field theory, where one has (3 = | and 5 = 3, this correlation length is equivalent to 
that one derived from the second moment of the correlation function, which in momentum 
space is given by G^q) - 1 



From this relation one infers 33 



f = (r + 3m 2 ) 



'1/2 



with nib obtained from the mean-field equation of state: 

H = r nib + ml . 



(4.16) 



(4.17) 



Equations ([4.16| ) and ([4.1 7| ) lead to the mean-field expressions of the correlation lengths 
£± and Ih in Eqs. (|4.1| ) and (|4.2|) , respectively. The correlation lengths are experimentally 
accessible by measuring the scattering structure factor. This allows one to determine the 
nonuniversal amplitudes £q and Iq in Eqs. ( |4.1|) and (|4.2| ). By measuring the coexistence 
curve Mb = ±B\t\@ the nonuniversal amplitude B is obtained. 
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However, the amplitude Iq m Eq. ( |4.2|) is fixed once £q and the nonuniversal amplitude 
B are known: 



•Aca,c) \S R 



where the first term in parantheses is nonuniversal - but it does not depend on the definition 
of the order parameter - and Q 2 and R x are universal amplitude ratios leading to 
{Q 2 /8R x Y h « 0.38 ind = 3. 



Using Eq. (|4.1|) and m& = y |r| one finds a relation between the order parameter m& in 
the present model and the order parameter M& defined in an actual experiment: 

m b = * (4.19) 
(v^)^ 5 

We note that in a magnetic system the bulk field if is directly accessible, while for liquids 
we use the correlation length In as a measure of the bulk field H according to Eq. ( |4.2| ). 
So far we have assumed that the coexistence curve close to the critical demixing point is 
symmetric with respect to the temperature axis. However, this might not be the case in a 
real liquid mixture, so that one would then have to use appropriate linear combinations of 
the thermodynamic variables. 



V. RESULTS FOR THE FILM GEOMETRY 



In the previous section [see Eq. (4.5) in Subsec. [VD|, within the Derjaguin approxima- 



tion, we have used results obtained for the film geometry in order to analyze the limiting 
case of the colloidal particles being very close to each other. Besides that the film geom- 



etry is of considerable interest in its own right as the generic case for simulations |65|j66 



and for experimental studies involving wetting films p|,|3~0"| or force microscopes with crossed 
cylinders of large radii of curvature. The Casimir forces in the film geometry have been 



considered for walls that do not break the symmetry of the order parameter 25] and also for 



symmetry breaking walls |27] . However, most studies focus on the temperature dependence 
only, while in the following we consider also the case H ^ 0. 



For the film geometry the same numerical methods can be used as described in Sec, pi. 
However, since the system is effectively one-dimensional there are additional methods avail- 
able. For instance, the first integral of the second-order differential equation (|3.1|) can be 
found analytically |2"7fl . The integration constant, which follows from fulfilling the boundary 
conditions is proportional to the force (see App. A in ref. [27]). For H = the second 
integration leads to elliptic integrals, while for H ^ it must be carried out numerically. 

First we obtain the scaling function K±(^j, j^j of the singular contribution to the force 

between two parallel plates at distance L for the various thermodynamic paths shown in 
Fig.|2j. The temperature dependence (solid path) of the force is discussed in refs. 

'see in particular Fig. 3 in ref. [j4q| ) and reveals a maximum of the force at G r 



27 and 



L/£_|_ = 1.94. As a new result the dependence of the force on the bulk field H (dashed path) 



is shown in Fig. 10 in terms of the scaling variable H = sgn(if) L/Ih- 
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It is remarkable that the force maximum is about ten times larger than the force at the 
critical point. Since this scaling function enters into the expression for the force between 
two spheres this peak structure carries over to that case, too. Indeed, Fig. ^ shows that the 
force maximum in the limit of two very close spheres is about five times higher than the 
force at the critical point. 

Figure shows the dependence of the force on the bulk field for T ^ T c . Above T c one 
finds a similar structure like at T c except that the maximum is less pronounced. (Here we 
refer to the maximum value of the modulus of the force. However, in order to facilitate an 
easy comparison with the results of ref. |J7[] in Fig. |TT| the sign of the force is chosen such 
that it reveals its attractive nature.) Slightly below T c the maximum is more pronounced 
although the form of the curve is still the same. But at still lower temperatures the character 
of the curve changes. Instead of a smooth increase upon increasing the bulk field towards 
H = there is a jump from a small value of the force to a high one. This is accompanied 
by a first-order transition of the corresponding order parameter profile. 

Figure shows the two coexisting profiles at the critical field of this first-order phase 
transition for various temperatures T < T c . The profiles, which are negative in the middle of 
the film, correspond to the weaker of the two respective forces. This means that the forces 
are weaker if the solvent component preferred by the surfaces is the one in which the bulk 
solution is poor. These first-order phase transitions correspond to capillary condensation in 
the slit which can occur for fixed temperature as function of H as shown in Figs. [IT] and |TJ 
or as function of temperature for fixed values of H. The loci of these capillary condensation 
phase transitions are shown in Fig. |TJ. The transition line ends in the capillary critical point 
(©_ = 7.17, H = —9.54) which is in quantitative accordance with the findings for the critical 



point shift [68 



In ref. JJ7| the solvation forces between two parallel plates confining a two-dimensional 
Ising spin system have been calculated. They exhibit the same features as those found in 
the present mean-field calculations valid for spatial dimensions d > 4. For temperatures 
far below T c upon increasing the field the force jumps from a small to a large absolute 
value and the varies linearly as function of the bulk field (compare Fig. [IT]). In the linear 
regime the solvation force is given approximately by K^/L 4 = 2m(T)H where the 
bulk spontaneous magnetization fh(T), i.e., the magnetization at H = is fn = \t\ 1 / 2 within 



mean-field approximation which indeed gives the slope of the linear parts in Fig. |TT|. We 



note that in d = 2 capillary condensation can be identified only as a quasi-phase-transition 
whereas in d — 4 there is indeed a bona fide first-order phase transition. Figure [14] shows 
the difference 5K_ between the two corresponding force values at the transition. Upon 
approaching the capillary critical point this difference vanishes. Along the coexistence line 
H(T) of capillary condensation the quantity 2m(T) H(T) provides a good account of the 
force difference (see Fig. 14). 

The findings in d = 2 and d = 4 agree to a large extent qualitatively. Quantitatively, 
however, the results in d = 2 and d = 4 differ significantly. For example the ratio of the 
maximum force and of the force at the critical point as function of the bulk field and for 
T = T c is K l l(0, H max )/i4(0, 0) « 100 in d = 2 and ^(0, H max )/i^(0, 0) « 10.1 in d = 4. 
A simple linear interpolation yields an estimated value of if+(0, H max )/iCj|(0, 0) ~ 55 for 
this ratio in d = 3. Thus the studies in ref. [^7| and the present ones can be used to obtain 
quantitative estimates for the actual behavior in d = 3. 
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VI. CONCLUDING REMARKS AND SUMMARY 



We have studied the critical adsorption profiles and the effective free energy of interaction 
for a pair of colloidal particles which are immersed in a binary liquid mixture near its 
critical demixing point. Our results pertain to the whole vicinity of the critical point of 
the binary liquid mixture, i.e., they are given as function of both the reduced temperature 



t = (T — T c )/T c and the field h conjugate to the order parameter (see Fig.y and Sec. pJ\). 

In an ensemble of colloidal particles dissolved in a near-critical solvent, the critical 
Casimir forces between them are expected to lead to flocculation of the particles. This 
holds even in cases where dispersion forces alone are not strong enough to produce such 
a phase transition (see Subsec. [l"VF| ). Indeed, such a flocculation of colloidal particles has 
been observed experimentally for various binary liquid mixtures In particular, these 

experiments exhibit an asymmetry in the shape of the observed flocculation diagrams in 
that flocculation occurs on that side of the phase diagram where the binary liquid mixture 
is poor in the component preferred by the colloids. Our results are consistent with this 
observation, since we indeed find that the critical Casimir forces are strongest in this region 
of the phase diagram (see Subsecs. [lV B| , |1V Cj , and [IV E| ). 



In the following we summarize the main results of this work: 

1. For two spherical particles of radius R at a distance D (see Fig. [Tj) we have investigated 
the effective interaction mediated by the solvent, such as a binary liquid mixture. This 
effective interaction depends on the thermodynamic state of the solvent, for which the 
phase diagram is sketched in Fig.|2]. 

2. The solvent is described by a field-theoretic model, provided by the Ginzburg-Landau 
Hamiltonian [Eq. (|2.1|) 1, with additional terms representing the solute spheres. We 
solve this model within mean-field theory as described in Sec. |II[ 

3. ( [Ill A| ) The profile at the critical point is obtained by a conformal mapping of the 
corresponding profiles between concentric spheres. ( |III B| ) Close to the surface of either 
sphere the short distance expansion of the profile is valid, which is carried out off the 
critical point. ( [Ill C| ) For the determination of the complete profiles we first discretize 
the space using the conformal lattice shown in Fig. |3] and then apply the method of 
steepest descent for the order parameter at the lattice points. ( |III D| ) From the order 
parameter profiles the force between the spheres is calculated by integrating the stress 
tensor. The choice of this integration path is indicated in Fig. 01. 



4. The results for the two-particle interaction are presented in Sec. [TV]. (|I V A|) The de- 
pendence of the force on the temperature at the critical composition is shown in Fig. 
It exhibits a maximum above T c which approaches that of the critical point for large 



distances. ( [IV B ) The dependence of the force on the composition at the critical tem- 



perature is shown in Fig.y. It exhibits a pronounced maximum whose position also 
approaches that of the critical point for large distances. The values of the maximum 
of the force determined for various distances are shown in Fig.|7|. For small distances 
these values differ significantly from the value at the critical point. ( |IVQ) In Fig. || the 
dependence of the force on the temperature for fixed compositions off the critical one is 
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shown. ( [IV D|) Using the knowledge of the force between parallel plates the Derjaguin 
approximation can be applied for small particle distances. ( [IV E| ) The small radius 
expansion can be applied for large particle distances regarding the spheres as pertur- 
bation of the bulk phase. This expansion yields nonperturbative results for d = 3. 
In leading order the small radius expansion leads to forces which exhibit a maximum 
only at the critical point. (|I V F|) In Fig.|9| a comparison of the Casimir forces and 



the dispersion forces at T c and in d = 3 is shown. For large distances the dispersion 
forces decay faster while for small distances the amplitudes are smaller than those 
for the Casimir forces. Moreover, the temperature and composition dependences of 
the Casimir forces allow one to distinguish them from the dispersion forces. ( IV G|) 



Relations between thermodynamic quantities as they are used in the present model 
and their experimental counterparts are discussed. 

5. In Sec.[V| results for the case of parallel plates are presented, which are used via the 
Derjaguin approximation for the limiting case of small particle distances. Moreover, 
the present mean-field results are in qualitative agreement with corresponding results 
for a two-dimensional Ising spin system obtained in ref. [67J. Figure |TU] describes the 



dependence of the force on the bulk field H at the critical temperature. In Fig. |TT| the 
dependence of the force on the bulk field is shown also for temperatures considerably 
below T c . As H is increased the force jumps from a small absolute value to a large 
one. This is accompanied by a first-order phase transition of the corresponding order 



parameter profiles as shown in Fig. [12]. These first-order phase transitions correspond 
to capillary condensation in the slit and their loci are displayed in Fig. [13|. The transi- 
tion line ends in the capillary critical point. In Fig. |TJ] the size of the jump of the force 
at capillary condensation (see Fig. [I]]) is plotted as function of the bulk field. Away 
from the capillary critical point this difference is well approximated by an estimate in 
which only the order parameter value for H = at the transition temperature and the 
corresponding bulk field enter. 
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FIG. 1. The geometry of two equal neighboring spheres. 
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C A,c C A 

FIG. 2. Schematic phase diagram of a binary liquid mixture consisting of A and B particles in 
terms of temperature T and concentration ca of A particles. The shaded area is the two phase region 
— ending in the critical point CP — which separates the A-rich and the B-rich phase. The relevant 
thermodynamic variables near the critical point are t and h. Three different thermodynamic paths 
are indicated as considered in the main text. 
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FIG. 4. Calculating the force via Eq. ( 3.25| ). R is the radius of a sphere Vs that contains one 
of the spherical particles of diameter 2R. The projection onto the />z-plane is parametrized with 
the angle (p. The interior of V$ is shifted by the infinitesimal amount a, which is equivalent to an 
infinitesimal change of the distance D between the two particles. 
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FIG. 3. The conformal lattice (bottom) for two spheres of diameter 2 R at distance D (see Fig. [TJ) 
as obtained from the lattice between concentric spheres of diameter 2R+ and 2R-, respectively. 
(For the purpose of clear visibility here the two lattices are shown on different scales.) The sphere 
with diameter 2 Rq = 2 y/~K+ R- is mapped onto the vertical midplane of the conformal lattice. 
We note that inside the sphere with diameter 2 Rq the spherical lattice lines are chosen equidistant 
and more dense than outside. R and D are given by Eqs. (|3.9|) and ( 3.10| ). Vs is introduced in 



Subsec. HID 
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FIG. 5. Force between two particles as function of temperature T > T c and for H = relative 
to the force at criticality. The force is maximal for + = max (A) as indicated by the dots. 
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FIG. 6. Force between two particles as function of the bulk field H at T = T c relative to the 
force at criticality (r, H) = (0, 0). The force is maximal for 5 = H max (A) as indicated by the dots. 
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10~ 1 10° 10 1 10 2 10 3 10 4 
A = D / R 



FIG. 7. Maximum force values as function of distance. For small distances the force in- 
creases as A -5 / 2 which is scaled out here. Accordingly, the forces are normalized by the value 
Ko = limA— >o(A 5 ' 2 if+(A, 0, 0)) at the critical point in the limit of vanishing distance. For large 
distances the forces decay as A _2/3 /^ _1 , i.e., as A~ 3 within mean-field theory so that A 5//2 X + 
decays as A" 1 / 2 . 0+ = and 3 = sgn(H)D/l H . 
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FIG. 8. Force between two particles as function of temperature for fixed compositions ca ^ ca, c 
which implies fixed values of nib* = rrn,D^/ u . Considering the dotted path in Fig.||, i.e., for 
m-fc = const [for a relation between mj, and ~ (ca—ca,c) see, c.f., Eq. ( [4.19| )1 both thermodynamic 
variables r and H and thus both scaling variables 0+ = -£>/£+ an d 3 = sgn(H)D/lH are varied 
simultaneously [see, c.f., Eq. 4.17]. The Casimir force increases by going further away from the 
critical composition. 
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FIG. 9. Comparison of dispersion forces and Casimir forces at the critical point 
(T = T c , ca = ca,c) as function of the distance D in units of the sphere radius R. The forces 
are shown for both a sphere in front of a planar wall (SW) and for two spherical particles (SS). 
For a detailed discussion see the main text. 




FIG. 10. Normalized scaling function of the force for a film with thickness L at T c , i.e., 
G+ = L/£ + = as function of the scaling variable S = sgn(il) L/Ih- 
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FIG. 11. Normalized scaling function K± of the force between parallel plates at distance L for 
different temperatures, i.e., different values of the scaling variable 0± = L/£± as function of the 
bulk field H. Note that the scaling variable E = sgn(iJ) L/ljj raised to the power 5, where 5 = 3 
within mean-field theory, is directly proportional to H revealing the almost linear dependence of 
the force on H over a wide range of values. For details see the main text. 
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FIG. 12. Pairs of density profiles coexisting at capillary condensation for various temperatures 
T < T c relative to the bulk order parameter m^ p at the capillary critical point. The lower curves 
are stable further away from bulk coexistence whereas for the upper ones capillary condensation 
has taken place. Above the critical point of capillary condensation (see Fig.|l3|) the profiles look 
similar to the upper ones. 
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FIG. 13. Phase diagram for capillary condensation as obtained within the present mean field 
approach. The transition line ends in the capillary critical point (6- = 7.17, S = —9.54). The 
straight line ending at (6- = 0, S = 0) is the line of first-order phase transitions in the bulk. 




FIG. 14. The difference SK_ between the force values of profiles coexisting at capillary conden- 
sation (solid line) and an approximate expression thereof (dashed line). The overall L dependence 
of the forces (L~ 4 ) is scaled out in k][, 5K^i, and the approximation. 
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